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Long linear polymers in strongly disordered media are well described by self-avoiding walks 
(SAWs) on percolation clusters and a lot can be learned about the statistics of these polymers 
by studying the length-distribution of SAWs on percolation clusters. This distribution encompasses 
to distinct averages, viz. the average over the conformations of the underlying cluster and the SAW- 
conformations. For the latter average, there are two basic options one being static and one being 
kinetic. It is well known for static averaging that if the disorder of the underlying medium is weak, 
this disorder is redundant in the sense the renormalization group, i.e., differences to the ordered 
case appear merely in non-universal quantities. Using dynamical field theory, we show that the 
same holds true for kinetic averaging. Our main focus, however, lies on strong disorder, i.e., the 
medium being close to the percolation point, where disorder is relevant. Employing a field theory 
for the nonlinear random resistor network in conjunction with a real- world interpretation of the 
corresponding Feynman diagrams, we calculate the scaling exponents for the shortest, the longest 
and the mean or average SAW to 2-loop order. In addition, we calculate to 2-loop order the entire 
family of multifractal exponents that governs the moments of the the statistical weights of the el- 
ementary constituents (bonds or sites of the underlying fractal cluster) contributing to the SAWs. 
Our RG analysis reveals that kinetic averaging leads to renormalizability whereas static averaging 
does not, and hence, we argue that the latter does not lead to a well-defined scaling limit. We 
discuss the possible implications of this finding for experiments and numerical simulations which 
have produced wide-spread results for the exponent of the average SAW. To corroborate our results, 
we also study the well-known Meir-Harris model for SAWs on percolation clusters. We demonstrate 
that the Meir-Harris model leads back up to 2-loop order to the renormalizable real world formula- 
tion with kinetic averaging if the replica limit is consistently performed at the first possible instant 
in the course of the calculation. 

PACS numbers: 64.60.ah, 61.41. -|-c, 64.60.al, 64.60.ae 
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I. INTRODUCTION 



Linear polymers in disordered media have been an im- 
portant topic for experimental and theoretical study for 
more than 20 years It is a well known fact that 
that the universal scaling properties of linear polymers 
in strongly disordered media are well described by the 
statistics of self-avoiding-walks (SAWs) on percolation 
clusters. Usually, the term SAW implicitly refers to a 
mean or average SAW for which an average is taken over 
the number of steps or intrinsic lengths of all self-avoiding 
walks with a specified Eukledian start-to-end distance or, 
respectively, over the start-to-end distances of all self- 
avoiding walks with a specified intrinsic length. How- 
ever, there are other SAWs that are equally interesting. 
Most notably, there is the shortest and longest SAW for 
which the intrinsic length for given terminal-separation 
is shortest and longest, respectively. On critical perco- 
lation clusters, all these SAWs are fractals, i.e., their 
masses (which are proportional to their intrinsic lengths, 
of course) as functions of the start-to-end distance scale 
with non-integer scaling exponents. Interestingly, how- 
ever, it turns out that the mean SAW is more than just 
a simple fractal - it is a multifractal. On critical perco- 



lation clusters, the statistical weights of the elementary 
constituents (bonds or sites) are non-trivial, and an en- 
tire family of multifractal scaling exponents is required 
to characterize the distribution of these weights through 
its moments. 

Although the last 2 decades have brought great ad- 
vancement in the understanding of linear polymers in 
disordered media, there are certain problems that have 
caused enduring controversy. For example it turned out 
that the famous Meir-Harris (MH) model Q which has 
long been standing as the only existing field theoretic 
model for studying average SAWs on percolation clus- 
ters has trouble with renormalizability [3, Another 
example is the puzzling fact that sophisticated numerical 
simulations by various groups have produced widespread 
results for the scaling exponent vsAW describing the mean 
length of the average SAW, see [l| • 

In this paper, we highlight that one has to be careful 
about the notion of average SAW if the disorder of the 
underlying medium is strong. Namely, there are essen- 
tially two qualitatively different ways of averaging over 
all SAWs between two connected sites for a given random 
configuration of a diluted lattice, one being static and the 
other being kinetic. It is well known that the exponent 
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i^SAW in a non-random medium is the same for static 
and kinetic averaging [^-Q , and that weak disorder of the 
medium in redundant in the sense of the renormahzation 
group (RG) if static averaging is used @ . Below, we em- 
ploy dynamical field theory to demonstrate that the same 
holds true for kinetic averaging. To discuss the effects of 
strong disorder, we present a renormalizable field theory 
for SAWs on percolation clusters based on the random 
resistor network (RRN). We use this theory to calculate 
the scaling exponents of the shortest, the longest and the 
average SAW as well as the entire family of multifractal 
exponents for SAWs on percolation clusters. This the- 
ory demonstrates that kinetic and static averaging may 
lead to very different results if the random medium is 
at the percolation point. Within the real-world inter- 
pretation of Feynman-diagrams, static averaging leads 
to non-rcnornializability, and hence we argue that the 
statistics of linear polymers in disordered media may not 
have anasymptotic scaling limit, when static averaging is 
used. Because the static average has been used in many 
simulations, this might explain why the numerical results 
for vsAW are so wide spread. In fact, recent simulations 
by Blavatska and Janke [13, [HI using kinetic averaging 
are in excellent agreement with our theoretical results. 
To further corroborate our findings and to shed some 
light on them from a different angle, we include a discus- 
sion of the MH model in the context of kinetic vs. static 
averaging. A brief account of the work presented here 
has been given previously in Ref. [T2j . 



II. OBSERVABLES AND AVERAGES 

The fundamental question addressed in this paper is 
that of the scaling behavior the length of a (shortest, 
longest or mean) SAW on a percolation cluster when av- 
eraged over all cluster conformation. To this end, one can 
either consider the scaling of the mean Euclidian distance 
Rl between the starting point and a random endpoint of 
an SAW of length L, 



(2.1) 



or one can study averages over the random length L(x, y) 
of a SAW (which is proportional to the number of 
monomers, the intrinsic length or the mass of the cor- 
responding polymer) between a pair of sites {x,y) as a 
function of their Euklidian distance |x — y|. The latter 
approach is more convenient from the standpoint of field 
theory [3, and we will take it here. 

Let x{x,y) be the pair connectedness indicator func- 
tion which is unity if x and y are connected in the 
given random configuration C and zero otherwise. Let 
Z{w; X, y, €) be the generating function of the SAWs 7 
with length ^(7) = L belonging to the bundle B{x, y; £) 
of all SAWs starting at y and ending at x both on cluster 



C This generating function can be written as 
Z(zi;;a:,2/,e:)=^e-'^^ZL(a;,y,e:) 



L 

E 



p(7)e 



(2.2) 



where ^(7) is a weight function that depends on the av- 
eraging procedure one uses and which, of course, has to 
satisfy J^-yPil) — 1- Special cases arep(7) = if 7™ 
is either the shortest or the longest SAW belonging to the 
bundle B{x,y;<t). The perhaps most basic averaging pro- 

(B) 

cedures are either static, i.e., all the Ng/^ 7's belonging 

to B{x,y;€) are weighted equally, ^(7) = l/A^g^iV' 
kinetic, i.e., a given 7 earns a factor 1/z contributing to 
^(7) at each ramification where z — 1 other SAWs from 
the bundle B{x, y; £) split off and which, in general, leads 
to different weights for different 7's. Note that in kinetic 
averaging, the probabilities satisfy the additivity prop- 
erty p(7i U 72) = p(7i) +^(72) when two SAWs 71 and 72 
are identified, e.g., in a coarse-graining procedure. To the 
contrary, this is not the case in static averaging because 
the number Ng'^2w changes when SAWs are identified, 
and this fact leads to trouble when coarse-graining proce- 
dures are applied spatial inhomogeneous fractals like the 
backbone of a percolation cluster as they are in renor- 
malized field theory. 

In terms of the generating function, the mean length 
of the SAWs is given by 



d [x{x,y)hiZ{w]x,y,€)\ 



dw 



[xix,y)] 



(2.3) 



where [■ ■ ■ denotes an average over the configurations 
C in which each bond is occupied with probability p. 
At criticality, one expects scaling behavior of the mean 
length M[x, y) of long polymer chains, 



M{x,y) = [L{x,y)\^ 



|x- y 



1/ I'SAW 



(2.4) 



For kinetic averaging, in particular, the critical value Wc 
of w is zero. 

It is well known that multifractality can arise when 
physical processes unfold on fractals such as critical per- 
colation clusters. Typical examples are electrical con- 
duction on RRNs [l3l - [l8l | and random resistor diode net- 
works [H, [lOl where the distribution of currents flowing 
through the bonds is multifractal, i.e., is characterized 
by an infinite set of critical exponents which are not re- 
lated in a simple linear or affine fashion. It turns out that 
the situation is similar for SAWs on percolation clusters, 
where the moments 



with Sf, the length of bond h and 

mb= ^ Xb{l)p{l) < 1 

',<^B(x,yX) 



(2.5) 



(2.6) 
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the statistical weight of bond b with Xbil) = 1 if 6 belongs 
to the SAW 7 and Xbil) = if it does not, probe distinct 
substructures of the underlying percolation cluster and 
hence the average 



\Y(x,y)L''°'Hx,y)] 



(2.7) 



leads to an infinite family of multifractal exponents v^""); 

MM(x,2/)^|x-y|i/'''°\ (2.8) 

One has the special cases v^^^ = 1/Dbf,, v^^^ = vsAW, and 
j^(oo) _ ijjj^^j^ — u where Dbb is the fractal dimension 
of the backbone, Dred is the fractal dimension of the 
red (simply connecting) bonds, and v is the percolation 
correlation length exponent. 



III. SAWS IN DISORDERED MEDIA AS A 
KINETIC PROCESS - THE EFFECTS OF WEAK 
DISORDER 

In this section we briefly discuss the effects that weak 
disorder of the underlying medium has on the statistics 
of SAWs. For static averaging over SAW conformations, 
it well known that this disorder is redundant in the sense 
of the RG, i.e., it affects only non-universal quantities 
(in annealed as well as in quenched disorder-averages) 
although a nai've application of the Harris criterion [9[ 
apparently signals its relevance. This redundancy can be 
shown by introducing non-correlated quenched disorder 
into the usual (/9^-field theory of a n-component order pa- 
rameter which, when the replica trick with m-fold repli- 
cation is used to facilitate averaging over the quenched 
disorder, generates the statistics of SAWs with static av- 
eraging via the limits 71 — > and m — > 0. Here, we are 
interested mainly in kinetic averaging, and hence it is of 
some interest to demonstrate that this redundancy also 
holds for this type of averaging. A field theory for ki- 
netically generated SAWs in ordered media has been for- 
mulated by Pcliti more than twenty years ago Q. Here, 
we introduce a dynamical field theoretic model for kineti- 
cally generated SAWs in disordered media, and we utilize 
this model to discuss the effects of weak disorder. To set 
the stage, we first introduce reaction diffusion processes 
that model SAWs in disordered media as a kinetic process 
and then we jump to the dynamical response functional 
for these processes. Additional background and some de- 
tails of its derivation using the creation-destruction op- 
erator formalism are given in Appendix |^ It should 
not go un-noted that this description of SAWs in disor- 
dered media as a kinetic process completely avoids using 
the replica trick as well as the zero-component-limit and 
hence all the potential difficulties associated with these 
limits. 

On a very basic and intuitive level, kinetic SAWs in 
disordered media can be described by a set of simple 



diffusion and reaction processes. As in the problem of 
Brownian walks, the walkers hop on a d-dimensional cu- 
bic lattice from a site to a neighboring one + This 
is described by a random process 



A{y) 4 A{v + 8) 



(3.1) 



where A{y) denotes a random walker at r, and A is the 
hopping rate. Self-avoidance is introduced by production 
of markers B{y) at r, and destruction of the walkers by 
interaction with a marker 



^(r) 4 A{v) + B{y) , 
A{v) + B{r) A B{v) , 



(3.2a) 
(3.2b) 



where a and j3 are reaction rates. In addition to these 
processes, there is quenched disorder modelled by static 
traps C(r) with a random distribution p{v), and this dis- 
order acts on the walkers via a destruction process with 
rate 7 



^(r) + C(r) A C(r) , 



(3.3) 



Due to the markers and the traps, walkers living at time 
t should avoid sites that they have visited at times t' < t 
as well as sites which are locations of traps. 

These diffusion and reaction processes can be con- 
densed into a field theoretic functional following the work 
of Pcliti. A series of steps which are sketched in the Ap- 
pendix leads to the dynamical response functional 



J 



-J 














■/ —00 



dt s 



dt ss 



(3.4) 



s(x, t) is a field that has its origin in the variable A and 
encodes the random position x of a walker at time t ac- 
cording to the above reactions. Its mean value is the 
probability density finding a walker at these coordinates. 
s(x, t) is the corresponding response field that creates a 
walker at position x at time t. The dynamical response 
functional is invariant under the duality transformation 
s(x, <) o s(x, — t). A is the usual kinetic coefficient. The 
dependence of J' on the disorder is hidden in the pa- 
rameter T and in the coupling constant g. The latter 
consists of two parts: a positive part stemming from the 
self-avoidance and an additional negative part stemming 
from the disorder. Note that the fact that disorder and 
self-avoidance generate the same type of coupling in the 
field-theoretic functional makes the Harris criterion in- 
applicable for the problem at hand. 

Analyzing the RG fiow, it turns out that the fixed 
point value of g which corresponds to the limit of asymp- 
totically large SAWs is independent of the nonuniversal 
value of g and hence of the disorder fluctuations. Be- 
cause we are mainly interested in this limit, we can as- 
sume that the sole remaining disorder-dependence rests 
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in r. However, we can eliminate r from the dynami- 
cal response functional by letting s{t) s(t)exp(— At) 
and s{t) — > s{t) exp{Xt). This implies that the disorder- 
averaged Greens functions - the probability densities for 
finding TV walkers at positions {x} at times {t} if they 
are created at positions {y} at times {t'}, respectively - 

/ N N \ 

Gw({x, t}, {y, t'}) = / II s(x„, t,) Y[ ~s{yp,t'p) \ , 

\a=l /3=1 / 

(3.5) 

where (• • • )('^) denotes the cumulants with respect to the 
Boltzmann weight exp(— ^J) and denotes disorder av- 
eraging, are of the form 

G^({x,0,{y,i'}) = Gn ({x, t}, {y, t'})o 

, N . 

xexp AE(ia-ta) , (3.6) 

^ Q = l ^ 

where the index indicates the Greens functions without 
any disorder. The differences (^q — t^) are proportional to 
the lengths of the corresponding SAWs, and hence, the 
exponential factors correspond to a change in the non- 
universal fugacities in the statistics of the SAWs. Thus, 
Eq. p.6|) reveals that all universal properties remain un- 
affected by disorder. This establishes that the universal 
properties of kinetically generated and averaged SAWs 
are independent of weak disorder as their statically aver- 
aged counterparts are 

IV. NONLINEAR RANDOM RESISTOR 
NETWORKS 

The RRN is a variant of the usual percolation prob- 
lem where occupied bonds are viewed as resistors and 
open bonds are viewed as insulators. Here, we con- 
sider a nonlinear generalization nRRN of the RRN for 
which it is well known that the total resistance between 
2 points becomes proportional to the length of the short- 
est and longest SAW between those 2 points, respec- 
tively, for specific limits of the nonlinearity. In previous 
work [3, mi, US [13 , we have shown that the Feynman 
diagrams for RRNs (including their nonlinear generaliza- 
tion) have a real- world interpretation, i.e., they can be 
considered as being resistor networks themselves. Based 
on this real-world interpretation, we here develop an in- 
tuitive and powerful field-theoretic method to calculate 
the mean length and multifractal moments of SAWs on 
percolation clusters. 

To be specific, we consider a d-dimensional lattice 
where each bond is randomly occupied with probabil- 
ity p by a conductor or empty with probability 1 — p. At 
each lattice side i there is a voltage Vi . The voltage drop 
at a bond {ij) between sites i and j obeys a generalized 
Ohm's law 



with a bond resistance P[ij). Equivalently, with s = l/r 
the current /jj- across the bond is given by 

h,^a^,,)\V,~V^-' {V,-V,), (4.2) 

where = p^^-, is the non-linear conductance of the 
bond. We restrict ourselves in the following to the case 
that all occupied bonds have identical elementary con- 
ductances a ~ , and the conductances of the unoccu- 
pied bonds are zero. The currents are conserved at each 
site and obey Kirchhoff's first law 

^4,+/, = 0, (4.3) 

3 

where /; is an external current li flowing into site i. If 
there are only two ports x and y, these currents are given 
by li = I^Si^y — Si^x) , where / is the current resulting from 
the voltage difference U = Vy~Vx between the two ports. 
The electrical power P dissipated in the network is given 
by the bilinear form 

= E'^fe) - ^»r^' = E/'fe) i^^^^r^' > (4-4) 

where the two last equalities follow from Ohm's law, 
Eqs. dHI) and dOl). Using Eq. (gj]), one has 

P = C// = i?,(x,2/)|/r+^ (4.5) 

in the two-port case, where Rr{x, y) denotes the total re- 
sistance of the network between the two ports. Blumcn- 
fcld et al. p6j have shown that the special values r = — oo, 
— 1, —0, -|-0, 1 and oo describe physically relevant geo- 
metric properties of the diluted lattice. In particular, it 
is easily demonstrated that for r — > ±0 the internal cur- 
rents at each ramification flow only in the direction of the 
highest [r — > -1-0) or the lowest (r — > —0) voltage gra- 
dient (electromotorical force) thereby mapping out the 
shortest or the longest SAW between the 2 terminals, 
respectively. As a consequence, the resistance Rr{x,y) 
is proportional to the length of shortest or longest SAW 
between x and y for r — -|-0 or r — > —0, respectively. 
We will use this fact for calculating the average length 
M{x,y) for these SAWs via calculating the total nonlin- 
ear resistance Rr{x, y) averaged subject to the condition, 
that the two ports are on the same cluster, 

Mr{x,y) = [x{x,y)Rr{x,y)\J[x{x,y)\^. (4.6) 

At criticality, the average total nonlinear resistance obeys 
the power law 

M,(a;,y)~|x-y|i/^- . (4.7) 

which will allow us to extract the SAW exponents for 
the shortest and the longest SAW simply by taking the 
appropriate limit with respect to r. 
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As far as the average SAW is concerned, the situa- 
tion is somewhat more subtle. Obviously, the average 
length of the average SAW lies in between the average 
length of the shortest and the longest SAW, which are, 
of course, very different. Since average SAW sits some- 
where in this discontinuity at r = 0, it is not known how 
to extract its average length from the nRRN by a limit 
taking. To overcome this problem, we developed the idea 
to study the average SAW by using the real- world inter- 
pretation [l^, [2^ [23, of Feynman diagrams which 
we will discuss in detail in the following section. For 
studying SAWs on percolation clusters, we extend the 
real-world interpretation originally developed for study- 
ing electrical transport on RRNs in that we put SAWs on 
Feynman diagrams. As we will explain in detail below, 
the task of calculating the average length of SAWs on 
percolation clusters then in essence reduces to calculat- 
ing the average length of SAWs on Feynman diagrams. 
For the average SAW in particular, this approach avoids 
the aforementioned problems associated with taking a 
limit in r 

The resistance Rr{x, y) can be obtained by solving the 
circuit equations (|4.2I) , (|4.1I) , and (|4.3p . The circuit equa- 
tions can be viewed as a consequence of the variation 
principle 



d 
dVi 



1 



s + 1 



P{{V}) - I(V, - Vy) 



0. 



(4.8) 



where the power P is expressed purely as a function of 
the set of all voltages {V}, see Eq. ()4.4|) . Obviously the 
network may contain closed loops. Suppose there is a 
complete set of independent currents {Z^'-*} circulating 
around these loops. Using Kirchhoff 's first law (|4.3|) and 
Eq. (|4.4p . one can express the electrical power P entirely 
as a function of the external current / and the set of loop 
currents {/(')}. Then, one readily obtains Kirchlioff's 
second law as a consequence of the second variational 
principle 



d 



dm 



P(/,{/('^})=0 



(4.9) 



This equation is used in the following to determine the 
loop currents as linear functions of the external current 
/. For r > 0, the variation principle (j4.9p has a unique 
solution which corresponds, of course, to the global min- 
imum of the power P. Without ambiguity, this solution 
leads to Rr^+Q{x,y) ^ L^in{x,y) where Lmin denotes 
the length of the shortest SAW (the chemical length). 
For r < 0; the situation is less straightforward. There 
exist in general several solutions which arc all local max- 
ima of the power P [2^ . Only if one selects the solution 
corresponding to the global maximum of P, one gets the 
length of the longest SAW via Rr^-o{x, y) ^ Lmax(a;, y). 
As mentioned above, the average length of the average 
SAWs has to lie somewhere in the interval between these 
two extremal values and for a correct interpretation of 
Ro{x,y), it seems natural to demand that it produces 
the length of the average SAW, Ro{x,y) ~ LsAw{x,y)- 



HARRIS MODEL 

A field theory for the non-linear random network was 
set up by Harris 27 1 in analogy to the field theory 
28 , Isof . The network is replicated 



of the linear case 

D fold: ^ = {V^ 



(1) 



V^ ^'). One considers 



the correlation function G{x,y;X) = {"i/ ^{x)'i> _^{y)) of 
^^X^x) cxp(iA • Vx) with complex currents iX ^ 0: 

G{x,y;X)= [z'^ | J] H ^^/"^ '^^P ( " TTT^^^^^^ 



j a=l 



+ iX-iVx-Vy)) . (5.1) 



Here P{{V}) = J2a PilV^"^}), and Z is the usual con- 
figuration dependent normalization. In order to be well 
defined, the integrations over replicated voltages are aug- 
mented with extra weight factors exp{iLuV^). Physically, 
these weights correspond to grounding all sites via unit 
capacitors. In this picture, lu with Imw > corresponds 
to the frequency of the voltages. 

Because the electrical power depends only on voltage 
differences, the integration over the mean voltage of each 
independent cluster of connected conductors leads to cur- 
rent conservation for this cluster in the limit a; — >■ 0. 
It follows that Xx = Xy ~ X ii the ports x and y are 
connected. However in the case that x and y are not 
connected, there arise factors ~ exp(cA^/iw), where c 
is some positive constant, which go to zero in the limit 
ioj — )• —0. Here, the condition A 7^ is essential. Then, 
as the result of this integration, the pair connectedness 
indicator function of the two ports xi^^ u) is automat- 
ically generated along with other factors which go to 
one in the limit D ^ Q. Following the work of Harris 
p7| , the integration over the voltage differences can now 
be done by the saddle-point approximation if we chose 
A(") = -il + ^(") with Y.a C'"^ = under the condi- 
tions 1 < p|/r+^ < D-^ and p\rr-'^^\ < 1 which in- 
deed can be satisfied simultaneously in the replica limit 
Z) — > 0. Note that the saddle-point equations are identi- 
cal with the variation principle stated in Eq. ()4.8p . Thus, 
the saddle-point is determined by the solution of the cir- 
cuit equations (|4.1[) . (|4.2p . and (|4.3|) . and according to 
Eqs. (|4.4p and (|4.5p . we obtain 

G{.x, y; X) = \xix, y) exp (^^Rr[x, y)) 

= [x{x,y)]^{l + ^Mr{x,y) + ...}, (5.2) 

where A,. = X^fLil^*'^'""'')'^^^ where we dropped 
a factor that goes to 1 in the limit D — >■ 0. Hence, 
G{x, y; A) is the cumulant generating function for the re- 
sistance Rr{x,y) between the connected ports x and y. 
Rr{x,y) is proportional to the elementary resistance p. 
Hence, limp_>o G(a;, y; A) = G{x,y;X ^- 0) = [x{x,y)]p 
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is the correlation (connectedness) function of the perco- 
lation problem. 

To safely exclude A = from the theory it is use- 
ful to resort to a lattice regularization of the voltage- 
integrals [2^ . One switches variables V to 9 = k and 

A = XVn ^ t^'^i'^S discrete values on a ZJ-dimensional 

torus, i.e. k and I arc chosen to be ZJ-dimensional inte- 
gers with -TV < /fc("),/(") < 7Vand/fc(") = fc(")mod(27V), 
= /t(")mod(27V). A is a redundant variable with ar- 
bitrary scaling behavior. 

After discretization there are {2N)^ — 1 independent 
state variables per lattice site, and one introduces the 
Potts-spins 

<^g{x) = (27V) e^p(*^ • ^ *a(^) = - c^^y" 

(5.3) 

subject to the condition J^g^e^^") ~ ^- essential as 
we already have remarked above that the limit I? — > is 
the first of all involved limits and, in particular, has to 
be taken before — > 00. 

The replication procedure leads to the effective Hamil- 
tonian 

-ffrcp = - In cxp ( —P 



1 



(5.4) 



which may be expanded in terms of or, equivalently. 



op 



J2 E^^'W*-A-(^)*A-(^') 

At^O 



E 



(5.5) 



<x.x'> a 



where is the (discrete) gradient in the replica space. 
Next the kernel K{X) is expanded in the limit of large 
conductance a (small resistance p) 



K{X) = A'o + A'lA, + K2AI 



(5.6) 



with K„ ^ p" . Therefore, in the limit cr — >■ 00 we 
have Ar(A) Kq, and we recover the (2A^)^-state Potts- 
model which describes percolation in the limit A* — >■ 0. 

By choosing the redundant variable A appropriately, 
one can show that all terms in Eq. (|5.6p with A"„>2 are 
irrelevant in the sense of the renormalization group. The 
Kn>2 merely lead to corrections to the leading scaling 
behavior that have been calculated for the linear RRN, 
r — >• 1, in Ref [sl]- They do not lead, however, to a 
family of crossover exponents as erroneously concluded 
in Refs. [10, HH. Bluntly stated, RRNs are multifractal 
but they are not multicritical. 



VI. FIELD THEORY 

To set up a field theoretic Hamiltonian T-L. one proceeds 
with the usual coarse graining step and replace the Potts- 



spins ^g{x) by the order-parameter field (/3(x, 0) defined 
on a d-dimensional spatial continuum. Constructing all 
possible relevant invariants of the symmetry transforma- 
tions of the model, performing a gradient expansion and 
discarding all irrelevant terms, one arrives at the Hamil- 
tonian 



1-L= dfix 



'■+1 , 9 3 



(6.1) 

Here r and w are the strongly relevant critical control pa- 
rameters and {-deY^^ := Y.ai-^ldO^'^^Y^^- For w = 
which corresponds to p = 0, one has full Potts-symmetry, 
that is invariance under all permutations of the sym- 
metric group S(^2N)°- The {2N)^ discrete states of the 
order-parameter field p{x,9) — "^^-^QCxpi^iX ■ 6) ^^{x) 
transform as the fundamental representation of this per- 
mutation group. This fact is crucial because it is the irre- 
ducibility of this representation that ensures that there is 
only one invariant second and only one invariant third or- 
der coupling with one unique relevant control parameter 
T and coupling constant g, respectively. 

Now we set up a diagrammatic expansion. Contribut- 
ing elements are the vertex with weight —g and the Gaus- 
sian propagator which reads 



l-%.o 
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'^A,0 



wAr{X) 



wAr{X) 



(6.2) 



in Fourier space. This inclusion-exclusion equation shows 
that the principal propagator decomposes into a propa- 
gator carrying A (conducting) and one not carrying A 
(insulating). This fact allows for a schematic decomposi- 
tion of the principal diagrams into sums of RRN-like con- 
ducting diagrams consisting of conducting and insulating 
propagators, see Fig. [1] Note that, like the momenta p, 
the (imaginary) currents A are conserved at each vertex 
of the conducting diagrams. For actual calculations, it 
is more practical to use continuous rather than discrete 
replica currents A. Once the decomposition is done, it 
is save to switch back to integrations over loop replica 
currents using 



A//V 



l[dX^" 



(6.3) 



where the limit 13 ^ is understood. 

The Feynman diagrams resulting from the decompo- 
sition have a simple and intuitive interpretation |2l| 
which is closely related to the link-node-blob picture of 
the back-bone of a percolation cluster. The diagrams 
may be viewed as being resistor networks themselves 
with conducting propagators corresponding to conduc- 
tors and insulating propagators to open bonds. In a 
Schwinger proper time parametrization where each con- 



7 




FIG. 1. Decomposition of self-energy diagrams to 2- loop or- 
der, and enumeration of the lines. 



ducting propagator i is written as 



1 



+ 'wAr{Xi) 



/ dsi expl- {pf+T+wAr{Xi)) Sj] 
Jo 

. ^^-^^ 

the Schwinger parameter Si corresponds to its resistance 
Pi — Si, and the rephca variables —iXi to currents flow- 
ing across this conductor. As mentioned above, this real- 
world interpretation is closely related to the link-node- 
blob picture. The conducting propagators correspond to 
the links, the vertices to nodes and the self-energy dia- 
grams to blobs which themselves can consist of a network 
of links with blobs. The backbone corresponds to the 
full Greens function, i.e., the propagator with all possi- 
ble self-energy insertions summed up by the Dyson equa- 
tion. When Fourier transformed back from momentum- 
space to configuration-space, the Schwinger parameter 
(the proper time) of a conducting propagator is propor- 
tional to the intrinsic length of a tortuous link generated 
by a diffusional motion, and thus the real- world interpre- 
tation naturally assigns to links their proper length. 

The conserved replica currents may be written as = 
Ai(A, {k}), where A is an external current applied at the 
external legs of the diagram (and subject to the Harris- 
conditions for application of the saddle point method) 
and {k} denotes the set of independent loop currents. 
The replica current dependent part of a diagram can be 
expressed in terms of its power P, 

exp ( — ^ SiAr(Ai) j =: exp[— ?«P(A, {k})] . (6.5) 



For the evaluation of the integrals over the independent 
loop currents we employ the saddle point method. Note 
that the saddle point equations constitute nothing more 
than the variation principle stated in Eq. (j4.9|) . Thus, 
solving the saddle point equations is equivalent to deter- 
mining the total resistance Rr{{si}) of a diagram, and 
the saddle point evaluation yields 



exp 



The Gaussian integration over all internal momenta 
is textbook matter. Thereafter, any self-energy diagram, 
see Fig. ([T]), or rather the mathematical expression stand- 
ing behind it is of the form 



/(p2. A) = /p(p2) - wIwip')AriX) -f 

l[ds, [l-wRri{s,})A,{X) 



^(p',{s.}), 
(6.7) 



where D{p^ , {si}) is the A-indcpendent part of the in- 
tegrand which is identical to the integrand of the corre- 
sponding diagram in the standard (/)^-theory. 

We use dimensional regularization and the renormal- 
ization scheme 

ifi (p = Z^^'^ifi , r — > f = Zt-t , 

w — ?> w = Z^^ZutW , 9^ 9^ = G^^ Z^^ ZuUfi'^ , 

(6.8) 

where e ~ 6 — (i, /i is an inverse external length scale, 
and Ge = (47r)-''/2r(l + e/2) is a factor that generi- 
cally emerges in the calculation of Feynman diagrams of 
a 0"^- theory. In the limit Z? — > 0, Z, Z,-, and Z„ are, 
the well-known percolation renormalizations calculated 
to three-loop order by de Alcantara Bonfim, Kirkham, 
and McKane [l^l . It remains to determine Z^, via calcu- 
lating the part of the self-energy diagrams that is propor- 
tional to w, see Eq. (|6.7p . Note that the renormalization 
factors have to fulfill several consistency checks for the 
field theory to be renormalizablc. Using the fact that 
the unrenormalized theory has to be independent of ^, 
one can set up in a routine fashion a Gell-Mann-Low RG 
equation 
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for the connected 
where 
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Gn, 



du 



dlnw 



dfj. 
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dfj. 
d\nZ 



dp, 



(6.10) 



3.11) 



are the corresponding Wilson functions. |o indicates that 
unrenormalized quantities are kept fixed while taking the 
derivatives with respect to p. Then, one can use standard 
methods to solve the RG equation at the infrared stable 
fixed point u*, determined by /3(u*) = 0. The Gell- 
Mann-Low function given to 2-loop order by 



/3(u) = pdf, 



(6.6) 
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(6.12) 
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leads to 



2 671 2 r^, 3, 



(6.13) 



for the fixed point. Augmenting the so-obtained solution 
with dimensional analysis, one gets the scaling form 

Gn (|x, wAr (jCj I ; T, M, /i^ 

(6.14) 

?7 = 7* and = (2 — k^,) ^ where = k(m*) efc. are the 
usual eritieal exponents for percolation. (f>r = (2 — Q^) 
is the resistance exponent. Choosing the flow parameter 
a,s £ — |x — x'|~^, Taylor-expanding the 2-point function 
G2 in powers of wAr{X) and comparing Eq. (j5.2p . one 
finally obtains 



Afr(x, x') ^ |x — x' 



'\'l>r-/l' 



.15) 



for the scaling behavior of the average total nonlinear 
resistance. 



VII. SAWS ON FEYNMAN DIAGRAMS 

In previous work, we have applied the real- world in- 
terpretation of Feynman diagrams to calculate the scal- 
ing properties of several physically relevant properties of 
percolation clusters: their average resistance when the 
bonds are linear resistors (r = 1), the fractal dimen- 
sions of the backbone (r — > —1), the minimal (chemical) 
length (r — > -fO), and the total length of the singly con- 
nected (red) bonds (r 00), as well as the multifractal 
moments of the current distribution [H, [HI, [H [H [11] . 
The key step in these studies was to determine the total 
linear or nonlinear resistance of the Feynman diagrams 
(or their multifractal moments in the study of multifrac- 
tality) as described above. In all these cases we verified 
that our theory was renormalizable. Furthermore, we 
have checked and verified our results were in conformity 
with results obtained by other methods as far as those 
exist. 

Now, we extend the real-world interpretation to study 
SAWs on percolation clusters. Instead of viewing them 
as networks on which electrical transport takes place, we 
view the Feynman diagrams now as media (or rather the 
backbones thereof) on which SAWs take place. In this 
picture, the conducting propagators correspond to links 
that are accessible to the walker and the insulating prop- 
agators are inaccessible. The Schwinger proper time pa- 
rameter Si of an accessible link corresponds to its internal 
curled length. The essential task is then to determine the 
(shortest, longest or average) total lengths 



(7.1) 



cf. Eq. p.Sp . of SAWs on the Feynman diagrams. The 
resulting mathematical form of the self-energy diagrams, 
in particular, is that of Eq. (|6.7p with Rr{{si}) replaced 
by L{{si}). From there on, after fixing the weights rrii 
of the propagators of the diagrams, the remaining cal- 
culation is once again textbook matter. Since L({si}) 
is a linear form of the Schwinger parameters, this cal- 
culation can be represented diagrammatically through 
self-energy diagrams with insertions into the conducting 
propagators. As indicated above, the length Lmin({si}) 
and imax({si}) of the shortest and longest SAW are pro- 
portional to Rr^+o{{si}) and Rr^-oi{si}), respectively. 
The length Lavc({si}) of the average SAW sits in the dis- 
continuity at r = and therefore can potentially provide 
helpful insights for its proper interpretation. 



A. The shortest SAW 

For calculating the average length of the shortest SAW 
on percolation a percolation cluster, we determine the 
total length of the shortest SAWs on Feynman diagrams. 
For a given self-energy diagram, that length is 



SAWs 



is SAWs 



(7.2) 



where the minimum is taken over all SAWs on conducting 
propagators connecting the external legs of that diagram. 
Details of the further steps leading from here to the expo- 
nent i^niin of the shortest SAW have been given in previ- 
ous publications [H, [H, [gj] , and we will not repeat them 
here. The upshot is that the diagrammatic expansion for 
the shortest SAW can be mapped onto that for dynam- 
ical percolation, at least to 2-loop order. This provides 
for an important consistency check for the real-world in- 
terpretation, and it provides also for a convenient way 
of calculating t'min by extracting it from the dynamical 
exponent z of dynamical percolation (s^ . The result is 



1 e 1231 45 /, „ 9 , „ 

2 + 24 + 2352 + 196 0'^2--ln3 



(7.3) 



B. The longest SAW 

In this section we calculate the scaling exponent i^max 
of the longest SAW on a percolation cluster. As detailed 
above, the length of the longest SAW between terminal 
points X and x' is proportional to the total nonlinear re- 
sistance between x and x' on that cluster in the limit 
r — >■ — 0. In the framework of the real- world interpre- 
tation, this means that we can calculate fc'max to 2-loop 
order via determining the total lengths 



-^max({Si}) 
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SAWs 



E ^ 

ieSAWs 



(7.4) 
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of the longest SAWs on the different self-energy diagrams 
depicted in Fig. [TJ Some Details of this calculation are 
presented in Appendix |B] It results in the renormaliza- 
tion factor 



u / 
^^Ie^\32e ' 128 
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192 
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(7.5) 

This result implies that the Wilson function is given 
by 
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Evaluating Ctu = 7 ^ 7«i the fixed point (|6.13p leads us 
then readily to our final result 
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for the inverse fractal dimension of the longest SAW. 



C. The average SAW 

As we have discussed above, we can apply the static 
or the kinetic rule to calculate the average total lenght 
L{{si}) of SAWs on a Feynman diagram. At one loop or- 
der, kinetic and static averaging lead to identical results. 
At two loop order, however, the situation changes, be- 
cause the 2 averaging procedures lead to different results 
for diagram H shown in Fig. [TJ Using the numeration of 
propagator- lines indicated in Fig. [1] the static rule leads 
to i^*^({si}) = (s2 -I- S4 -I- S5)/3 + 2(.si + 53)73, whereas 

the kinetic rule gives i^*"^({si}) = (si + S3 + S5)/2 + 
(s2 + 54)74. Having these 2 expressions for the aver- 
aged length, it is easy to understand that the static rule 
does not lead to a renormalizable theory. It is a basic 
fact of renormalization group theory that non-primitive 
divergencies arising from sub-integrations of a 1-loop in- 
sertion must be cancelled through the counter-terms in- 
troduced by the renormalization of this 1-loop insertion. 

However, the weights of i^*''({si}) are not compatible 
with the weights arising in the corresponding 1-loop di- 
agram with counter-term insertion: crunching the inser- 
tion to a point (corresponding to 52 + 54 — 0) leads 

to i^*^({si}) 55/3 + 2(51 + 53)/3 which is differ- 
ent from the total length of the 1-loop self-energy dia- 
gram with a point insertion. For kinetic averaging, how- 
ever, the additivity property mentioned in Sec. HIl comes 
into play, and crunching the insertion to a point gives 



L^'™''({5i}) — >■ (51 -1-52-1- 53)72 which is equal to the to- 



tal length of the 1-loop self-energy diagram with a point 
insertion. Hence, the kinetic rule produces non-primitive 
divergencies that are cancelled by the counter-terms from 
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TABLE I. Numerical values for various dimensions of the 
SAW exponents resulting from rational approximation. 



the 1-loop renormalization but the static rule does not. 
Thus, we have to reject the static rule on grounds of 
renormalizability, and we will use the kinetic rule in the 
following. 

The remaining steps in calculating the scaling expo- 
nent for the average SAW proceed as outlined above. For 
details, we refer to Appendix [Cl where the formulae for 
the multifractal moments reduce to those for the average 
SAW when we set the multifractal index a equal to 1. 
We obtain the renormalization factor 
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for the parameter w. Having Z^, it is straightforward to 
extract the SAW exponent vsASV as described above. We 
obtain the e-expansion 
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For comparison to experimental or numerical data, it 
is useful to improve the accuracy of our e-expansion by 
implementing rigorously known features. To this end, we 
craft rational approximations for z^min, ^^SAW, and t'max 
by adding 8th-order terms in e with coefficients chosen 
such that the exponents match the rigorously known re- 
sult J/. . . = 1 in d = 1 . Table U] compiles numerical val- 
ues resulting from this approximation for various dimen- 
sions. Figure [2] visualizes our £-expansions and rational 
approximations as functions of d. Note that our rational 
approximation for i^sAW agrees very well with the avail- 
able numerical estimates for this exponent which arc also 
shown in Fig. [2] 



D. Multifractality 

As mentioned above, the fascinating phenomenon of 
multifractality has been found in the past in situations 
where transport processes like electrical conduction take 
place on critical percolation clusters. It is reasonable 
to expect that multifractality also occurs in the con- 
text of SAWs on percolation clusters, and, indeed, it 
does [13, [13 • The length of the average SAW that we 
just computed corresponds to the first of the multifractal 
moments defined in Eq. (|2.5p . Now we allow the power a 
with which statistical weights of SAWs enter in Eq. (|2.5p 
to be arbitrary positive numbers. Doing so, we can influ- 
ence the way how SAWs contribute to the average over 
their bundle and thereby, loosely speaking, map out the 
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FIG. 2. The e-expansions of the exponents Umin (blue) [3^ . 
I'SAw (red), and Vnimx (green). The possible extrapolations 
at low dimensions d are shown by broken lines. The points 
denote compiled numerical results for vsaw P, [Ul ■ 



different fractal substructures of the mean SAW. This 
approach is guided by earlier work on RRNs, where mul- 
tifractality manifests itself in the moments of the current 
distribution. 

Using the real-world interpretation, we proceed in es- 
sentially the same way we did in Sec. IVIICl with the only 
difference that we now determine all moments of the sta- 
tistical weights of SAWs on the self-energy diagrams, i.e., 
we now keep the a in Eq. (|2.5p instead of restricting our- 
selves to a = 1. For details of this calculation, we refer 
the reader to Appendix [C] Here, we would like to point 
out, however, that this calculation further underscores 
the imperative of kinetic averaging because it leads to 
a renormalized theory even if all moments are included 
whereas static averaging does not. Our calculation pro- 
duces the renormalization constants 



1 



1 



47 
12 



1 \u 
2" / e " 
'2 1\ 1 
~ 2/ 4" 

This result implies the Wilson function 
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Evaluating C'") = 7 ~ t'"-* at the fixed point (|6.13p . and 
using — 1/(2 — (^i"') leads us then readily to our final 
result 
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for the family of multifractal scaling exponents defined by 
Eq. (|2.8p . As it should, this result for general a reduces in 
the special case a = 1 to our result for vsAW given above, 
and is perfectly consistent with the known results for the 
backbone and red bonds dimensions. This can easily be 
checked by setting a equal to and letting a — )■ oo, 
respectively. 



Note that Blavatska and Janke [lO|, |ll| have devised a 
Pade-type approximation of our e-expansion results for 
the multifractal exponents which comprise vsAW- This 
approximation agrees very nicely with their numerical 
results. 



VIII. MEIR-HARRIS MODEL 

In this section we discuss in some detail the RG of the 
Meir-Harris model for the average SAW on percolation 
clusters. Our motivation to do so is two-fold. First, we 
think that it is of some interest to shed light on the prob- 
lem at hand from a different angle, in particular, because 
we have no rigorous justification for our real-world inter- 
pretation based approach in the form of a mathematical 
proof. We will see below, that the MH model when renor- 
malized properly produces to 2-loop order the same result 
for the multifractal exponents as the real- world interpre- 
tation and hence provides a strong positive consistency 
check for the latter. Second, the RG of the MH model is 
very intricate and not properly understood even though 
the model has existed for more that 20 years now. A 
recent 2-loop calculation Q struggled with this intricacy 
and produced incorrect results. 

It is well known that the statistical properties of SAWs 
can be calculated from the m-component spin model with 
0(m)-symmctry in the limit m — > 0. To treat dilution, 
Meir and Harris start from the n-replicated version 
of the model. They introduce tensor fields \I'fe(x) = 

{"^k-lu'.'.'^aki^)}^ ^ < k < n, conjugate to the prod- 
uct of the replicated spin-components where the vector- 
indices ii are running from 1 to m and the replica indices 
ai ~ 1, . . . , n are arranged such that ai < • • • < a*;. Us- 
ing the Hubbard-Stratonovich transformation and pass- 
ing to the continuum limit, they obtain the effective 
Hamiltonian 



'J I, 



(8.1) 



Here, VE''^ is a symbolic notation for the sum over prod- 
ucts of three ^Pfe fields. Only those cubic terms are al- 
lowed for which all pairs (z, a) appear exactly twice. Di- 
agrammatically, this rule can be represented as shown 
in Fig. [3] for the ^^^3^2 couphng. Each of the SAW- 
representing replicons (thin lines) carries the indices of 
the corresponding field ^P. No two pairs of indices enter- 
ing an interaction vertex through a given inbound leg are 
permitted to exit the vertex through the same outbound 
leg. Furthermore, the SAW-limit m — > for the i-indices 
implies that diagrams in which some pairs of indices flow 
in closed loops produce vanishing contributions. Over 
all, any replicons flowing through an external line into a 
diagram must flow out off the diagram through another 
external line without making any internal loop. There- 
fore, the basic task is to count the different distributions 
of these self-avoiding replicons under the condition that 
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FIG. 3. Replicons flowing through a vertex. 

each line of the diagram bears at least one replicon. To 
circumvent the latter condition, it is useful to split each 
internal fat line (propagator with replica index fc > 0) 
into a difference of a conducting (fc > 0) and an insu- 
lating (fc = 0) line. This step leads to diagrams that 
can be drawn in the same way as those for the RRN, see 
Fig. [TJ After this decomposition, the next step is draw 
all possible self-avoiding replicons on the conducting di- 
agrams where, of course, replicons can flow only through 
conducting propagators. Then, one has to sum over all 
replicon distributions, i.e., all possible arrangements of 
internal replica indices by given external ones. Using el- 
ementary combinatorics, one finds that this summation 
for a diagram with N external legs produces a factor 

Zi{h,}) = lli^SAwi^,J)P ■ n (rfc' ^ ^ (8-2) 

where fcy = kji is the number of replicons entering at leg i 
(i = I, ■ ■ ■ , N) and exiting at leg j, and fc^ = fcy is the 
total number of replicons entering at leg i (i ~ 1, ■ ■ ■ , N). 
NsAwihj) is the number of different SAWs which can 
be drawn between the pair («, j) of external legs, ({^f^',}) 
is the multinomial coefficient fc/!/(fc/i! • • • fc;jv!)- Formula 
(|8.2[) reduces to iVg^vF ^'^^ self-energy diagrams with fc 
replicons and Nsaw different SAWs between the two 
legs. We parameterize the temperature-like control pa- 
rameters Ts by 

rs=f:(%i-r + sf2t}l:^^^+0{s'), (8.3) 

/=0 ^ ^ 1=1 

where t = vq- This parametrization facilitates the sum- 
mation over the replicon distributions after rj-insertions 
in the self-energy diagrams, as well as the limit fc < n — > 
0. The iij-part of a rs-inscrtion into the internal line p of 
an self-energy diagram, and summation over all distribu- 
tions of fc replicons leads to a factor 

Z(fc,/;p).,.iV|,^Q)(^^^)'.,, (8.4) 

where Nsaw{p) is the number of SAWs going through 
the line p. Hence, the ^/-insertions "measure" the Z-th 
power of the fraction of all SAWs drawn between the two 
external legs and going through the line p. In this sense 
the vi measure multifractal moments of the diagram using 
the static rule. 

However, the MH Hamiltonian (|8.ip is not multi- 
plicatively renormalizable as it stands. The order pa- 
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FIG. 4. 1-loop counter-terms. 




FIG. 5. 2-loop self-energy counter-terms. 



rameter fields belong to different irreducible tensor- 
representations of the direct product of the replica- 
permutation group Sn and the rotation group SO{m) 
for different fc. Hence, the fields ^fc need fc-dependcnt 
renormalization factors, and the model is critical at dif- 
ferent values (which are superficially set to zero in 
dimensional regularization) of the temperature-like con- 
trol parameters. Therefore, the model is highly multi- 
critical. For an earlier critique concerning this point see 
Le Doussal and Machta Furthermore, one needs in- 
dependent coupling constants gk.i,m for each product of 
three as opposed to a single coupling constant 

g because it is not possible to construct from the n-fold 
replicated m-vector model a higher simple symmetry- 
group where the order parameters for all fc belong to 
one and the same irreducible representation unlike in the 
case of the n-fold replicated m-state Potts-model leading 
to the m"-state Potts-model where such a construction 
is possible and commonly applied. The latter model, 
relevant for the dilute Ising model (to = 2) and the ran- 
dom resistor network (to — >■ 0) [s^l, therefore needs only 
one "scalar" coupling constant g and a unic (but to"- 
dependent) renormalization factor for all fields, and it is 
possible to apply the replica limit at the very end. It is 
not clear, however, for the MH model at which stage of 
its perturbation theory the replica limit should be taken. 
There has been hope that if the appropriate stage to 
take the replica limit can be identified the renormaliz- 
ability of the MH model in the form of a conventional 
multiplicative renormalization can be restored [38| . In 
the following, we will embark on a quest to identify the 
proper "timing" for the replica limit. 

One of the most basic facts of RG theory states that 
non-primitive divergences arising at a given loop-order 
in superficially divergent sub-diagrams must be canceled 
by the counter-terms resulting from lower loop-orders. 
The perhaps most direct route to understand this fact 
is provided by the iterative approach to constructing 
counter-terms invented at the dawn of RG theory by Bo- 
goliubov, Parasyuk, Hepp, and Zimmermann (BPHZ), 
see, e.g., Ref. [23 • We will use a BPHZ-like construction 
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FIG. 6. 2-loop vertex counter-terms. 



of counter-terms in a massless 2-loop calculation using 
t'Hoofts minimal dimensional renormalization (40| . For 
practical purposes, it is useful to split up the calculation 
into a part that determines the counter-terms of frames, 
i.e., those parts of Feynman diagrams that stand only for 
momentum integrations, and a part that determines dec- 
orations, i.e., symmetry factors, coupling constants and 
all other parameters that multiply the frames. There are 
7 frame-counter-terms to 2-loop order, sec Figs.|3]|ni Our 
calculation produces 



Sl^-^q\ VI = ^ 



S2-1 



3£ 

3e2 



q 

lie 



V2-1 



' 2e2 V 4 



6£2 V 12 



V2-3 



2e 



(8.5) 
(8.6) 
(8.7) 
(8.8) 
(8.9) 
(8.10) 



for these counter-terms. To calculate the corresponding 
decorations using formulas (|8.2p and (|8.4p . we determine 
all possible self-avoiding replicons on the conducting 1- 
loop diagrams shown in Figs. [1] and [7] as described above. 
We obtain for the 1-loop self-energy counter-term 



Selfl = ^(2'^- _ 2) • SI = ;^(2 - 2'')q^ 
2 ^ ' 6e ' 



5.11) 



where u = G^/i ^ ■ The Wi-insertions into these dia- 
grams result in the counter-terms 



Inslz- (■^^)wi = -g2^2^.-i_2 



"{2-2")ft)«, (M2) 



and the 1-loop vertex counter-term with A:, Z, and m repli- 






FIG. 7. Decomposition of the 1-Ioop vertex diagram. 

cons in the external legs is 

Vertl • g = {2^k+i+ra)l^ _ 3)^^^^^^ . yi 

= -(3-2('^-+'+™)/2)iV,,„„5, (8.13) 



where we have used the notation 

k\l\m\ 



2 • 



.14) 



f k+l—rn \ I { k+m—l \ \ ( m-\-l — k \ \ 

A 2 )-y 2 )-y 2 )■ 

For the MH model, we use the renormalization scheme 



1'^* = Z'/^4', (8.15) 
VI ^ VI = Z-^Zivi , (8.16) 
g^g^Z-^/^ZgCj. (8.17) 



To 1-loop order, the above counter terms are related to 
the renormalization factors introduced by this scheme via 



Z = 1 
Zi = 1 
Z, = l 



Selfl -\ 
Insl( H 
Vertl 



(8.18) 
(8.19) 
(8.20) 



In the replica limit n — > (which implies vanishing exter- 
nal replicon numbers fc, . . .) we retrieve the well-known 
percolation renormalization factors to 1-loop order. In 
particular, we retrieve 



1 



2u 

e 



^.21) 



Note that the renormalization of r follows from the iden- 
tity Zr ^ Zoo- 

Now we turn to the 2-loop counter-terms using the 
same general approach as for the 1-loop part of the cal- 
culation. From the diagrams shown in the middle part 
of Fig. [TJ we obtain the counter-term 



Self2-1 



2k 



6) • S2-1 



" fl_£V(4fe-5-2'= 
6e2V 3/ ^ 



6) - q. (8.22) 
The diagrams shown in the lower part of Fig. [T] yield 



a* 

Self2-2 = y (/i''^ - 3 • 2'' + 3) • S2-2 



u 



36e2V 12 



5.23) 



Note the placeholder h appearing in this formula. This 
placeholder reflects the fact that there are apparently 
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FIG. 8. Insertions into the 2-loop self-energy diagram whose 
decomposition is shown in the middle part of Fig. [T] 



^^^^^ 



replica limit as mentioned above. When we let n — >■ in 
the 1-loop superficially diverging self-energy subdiagram 
in diagram H before taking the summation over the repli- 
con distribution, we obtain a = 1/2, b = 1/4, c = 1/2 
and h ~ 2. Otherwise, we get a = 2/3, h = 1/3, c = 1/3 
and /i = 3. Wc will analyze the correct "timing" of the 
replica limit in more detail further below. 

To calculate the renormalization factors for the fields 
to 2-loop order, we collect our various diagrammatic re- 
sults, 



FIG. 9. Insertions into the 2-loop self-energy diagram whose 
decomposition is shown in the lower part of Fig. [T] 



Z =1 + Sclfl + Sclf2-1 + Sclf2-2 + . 
Zi = l + Insl( + Ins2-1; + Ins2-2i + 



(8.26) 
(8.27) 



two possible choices for taking the replica limit, and the 
result we obtain for diagram H of Fig. [1] depends on this 
choice. We can let n — )- in the superficially diverging 
1-loop self-energy subdiagram in H either before or after 
taking the summation over the replicon distribution. In 
the first case, h = 2 whereas ft, = 3 in the second. We 
will return to the issue of these choices further below. 

Next, we consider Ui>o-insertions into the 2-loop self- 
energy diagrams. From the diagrams shown in the middle 
part of Fig. [1] wc obtain 



Ins2-1, ■ \ ^]vi 



- ^ (4''- (1/2)' - 4 • 2'= • (1/2)' + 2) (^^ VI ■ V2-3 
2/(4*^' • (1/2)' - 3 •2'^' • (1/2)' - 2^^ + 3) Q^w; ■Y2-1 



2k-l r,2+fe-i 



- 2^*^"' - 2 



+ (22^-"' - 3 • 2*^-' - 2" + 3) (l - I) } • [Jj VI . (8.24) 

Insertions into the 2-loop self-energy diagrams shown in 
the lower part of Fig. [T] produce 



Ins2-2/ • yj^jvi 
= -g'^{h^ • a' - 2 • 2^= • (1/2)' - 2^ + 2) {^^vi ■ V2- 



.g*(ft'^'-&'-2-2''^(l/2)' + l) 



VI ■ V2-2 



- g^l^h" -2-2" ■ (1/2)' + 1) ( \ vi ■ V2-1 



J{-(/,''^a'-2i+^-'-2'^- + 2)i(l-I|) 



(K' -U -2 



l)-fl-- 
''2 V 4 



)} 



^.25) 



were a, &, and c are further /-independent placeholders 
stemming from the two apparent choices for taking the 



and take the limit fc 



Z = l + - 
6e 



37 

11 - — £ 



0. In this limit wc obtain 

2 



u 



12 /36e2 



Z, = 1 



1-2 



47 \ / 29 , 



(8.28) 



— e a 
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2 J V3 36 / 



u 

4^ 



5.29) 



As usual, these renormalization factors, as well as their 
products have the form of a Laurent series, Z = 1 + 

Our ultimate goal is to determine the inverse multi- 
fractal dimensions 



.(0 



= (2 - 



(8.30) 



of the MH model. Thus, we need to extract from the 
above renormalizations the Wilson function 



Ki^~f3{u)^hi{Z-^Zi] 
ou ' 



(8.31) 



where /3(u) = — eu + /3^°''(m) is the GcU-Mann-Low func- 
tion given to 2-loop order in Eq. (|6.12p . It follows from 
Eq. ((83T11 that 



^,=uUz-'Zif'-^-P^'\u)Uz-'Z,f^ 
au^ ' e ou ' 

1 d 



2e"Stx 



2{Z-^Zif^ -{{Z-^Zi 



(1) 



(8.32) 



has to be free of e-poles. Hence, we obtain the t'Hooft- 
identity ^ 



d_ 

du 



= /?(°)(.)£(z-zO 



2{Z-'Zif^ -{{Z-'Zif^ 



(1) 



5.33) 



Inserting our 2-loop results into this identity, we find the 
condition 



2a' - 65' + c' = 3 • 2"' - 6 • 4" 



.34) 
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This condition has the unique solution a = c = and 
b = 4^^. Thus, to make the Meir-Harris model renor- 
malizable, one necessarily has to take the replica limit 
n — >■ in the superficial divergent subdiagram (SDS) 
appearing in diagram H before one sums over the repli- 
con distributions of H with v;-insertion. Remarkably, the 
rcnormalization factor Zi obtained this way is identical 
to 2-loop order to the rcnormalization factor Z^j with 
a = I resulting from the Harris model in conjunction 
with the real-world interpretation provided that kinetic 
averaging is used. Consequentially, the same holds true 
for the multifractal exponents produced by the 

two approaches. We rate this as a strong indication for 
the validity of the real-world interpretation with kinetic 
averaging. 



IX. CONCLUDING REMARKS 

In summary, we have shown that weak disorder in the 
SAW-problem is redundant in the sense of the RG for 
kinetic averaging as it is for static averaging. We have 
derived the scaling exponents of SAWs in strongly disor- 
dered media by field-theoretic methods to second order 
in the dimensional expansion below six dimensions. We 
have shown in the real- world interpretation of the corre- 
sponding diagrams that in contrast to a static averaging 
over the SAWs only kinetic averaging lead to a renor- 
malizable theory. The different behavior of these two 
averaging procedures under renormalization is expected 
to have important physical consequences for the statis- 
tics of polymers in real disordered media. We argue that 
a statistics of polymers based on static averaging has no 
asymptotic scaling limit. Based on our findings, we do 
not expect experiments and numerical simulations using 
static averaging to produce clear scaling behavior. In 
fact, we think that the wide-spreading of simulation re- 
sults for the SAW exponent in strongly disordered media 
is linked to static averaging. 

Closing, we would like to supplement our firm but con- 
ceptually somewhat involved field-theoretic argument for 
the imperative of kinetic averaging by a simple hand- 
waving argument based on the link-node-blob model of 
percolation clusters. In this model, the backbone con- 
necting two terminal points of a percolation cluster, 
which is generically very inhomogeneous and asymmetric, 
can be envisaged as two nodes linked by a tortuous ribbon 
that contains blobs. A blob itself is constructed from at 
least two links joined at two nodes which may again con- 
tain blobs. Let us for simplicity consider an asymmetric 
blob as sketched in Fig. [10] that features two links be- 
tween two nodes, one with and the other without a blob. 
Note that this cluster resembles the ominous diagram H 
of Fig. [T] Assume that the internal blob has many rami- 
fications of short links in it. Hence, say A''' = — 1 ^ 1 
different SAW configurations are possible on the upper 
link. With static averaging the upper link acquires a 
much larger weight {N — 1)/N then the other (lower) 




FIG. 10. A blob in a blob of a percolation cluster in the 
link-node-blob picture. 

one (weight 1/A^) even if it may be much shorter than 
the link without the blob. Then, the statistics of the 
mean length is dominated by the short upper link with its 
many different SAWs induced by the blob. However, the 
weights change drastically upon coarse graining. Suppose 
we have some coarse graining procedure that culminates 
in condensing the "microscopic" blob into a single bond. 
After that, both links have the same weight. However, 
the lower one, since it is longer, now dominates the statis- 
tics. This demonstrates the instability of the weights of 
static averaging under real space renormalization as the 
group generated by repeated coarse graining. In contrast, 
kinetic averaging does assign the same weight to both 
links independent of the blob. Thus, kinetic averaging is 
stable under coarse graining even in a strongly asymmet- 
ric inhomogeneous disordered medium like the backbone 
of a percolation cluster. All in all, the behavior of the 
links- nodes-blobs cluster of Fig. [TOl under coarse graining 
resembles in a nut shell the issues we encountered in our 
discussion of diagram H of our field theory and thereby 
corroborates the imperative of kinetic averaging on an 
intuitive level. Note that real-space RG approaches as 
employed in Refs. [ll,^ generically use symmetric config- 
urations, and hence static and kinetic averaging lead to 
equal results in these approaches as they do in ordered 
media like regular lattices which are trivially homoge- 
neous and symmetric. 
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Appendix A: Dynamical response functional for 
SAWs in random matter 

In this appendix we generalize Peliti's derivation of 
a field theory for the statistics kinetically generated 
SAWs An excellent review of the method can 

be find in To have a starting point, let us re- 

visit the diffusion and reaction processes introduced in 
Sec. [ml First, the rules ^J^, ([S^) and for 
these processes are reformulated in terms of a mas- 
ter equation that describes the time dependence of the 
probability P({n, m}, i) for a given configuration of site 
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occupation numbers {n} = {. . . ,ni, . . .) and {m} = 
(. . . ,mi, . . .) of the walkers A and the markers re- 
spectively. Then, the master equation is transformed 
in the "second quantisation" formalism developed by 
Doi 43[, Rose |4J], and Grassberger and Sclieunert 
as follows. The configuration probability is encoded in 
the state vector \P{t)) I]{,i,m} ^(i*^' "^}' K*^' ™}) 
in a bosonic Fock space spanned by the basis \{n,m\). 
These vectors as well as the stochastic processes in 
the master equation are expressed through the ac- 
tion of bosonic creation and annihilation operators 
{«■*", &"*"} and {a, &}, respectively, which are defined 
via a+ |. . . , rii, . . . , {m}) = \. . . ,n^ + I, . . . ,{m}) and 
a^\. . . ,ni, . . . , {to}) = rii |. . . ,nj - 1, . . . , {m}) , etc. |0) 
is the vacuum state without walkers and markers. Subse- 
quently, the master equation can be written in the form 



-\P(t))=H\P{t)) , 



(Al) 



with an appropriate non-Hermitean pseudo-Hamilton op- 
erator 



H^xY,{a+- a+)a, + ^ {a(6+ - l)c 

<ij> i 

+ /3(1 - a+)a,h+b, + 7(1 - ci+)a,p.^ . (A2) 

Here < ij > denotes a pair of neighboring sites. The for- 
mal solution of the master equation (jAl[) reads \P{t)) = 
exTp{tH) |P(0)). Suppose we wish to calculate the 1- 
walker probability Pi{i,t) to find the walker at site i at 
time t if he starts from site at time not taking into ac- 
count the resulting distribution of the markers at t. Since 
a walker can be destroyed but not spontaneously gener- 
ated over the course of time, we have Pi{i,t) = {ni){t) 
where Ui = afui is the number-operator of walkers at 
site i and where (■ • • ) denotes an average whose precise 
definition will become clear shortly. To compute such a 
statistical average, it is useful to introduce the projec- 
tion state ('I = (0| niexp(aj; -I- bi). Using the identities 



and H \0) ~ one easily finds 



a+|0) = (.|a.c*^(a+-l)|0) . (A3) 



Changing from the occupation-number basis to a 
Bargmann-Fock space representation with the coherent 
states (the eigenstates of the annihilation operators) as 
the basis, and following standard procedures [1, IH, E^] 
the expectation value (|A3|) can be expressed as a path 
integral 



Pi{i,t) = J 23[a, a, S, fe]a,;(t)ao(0) exp(— 5[{a, a, 6, 6}]) 
=: (a,(t)ao(0)) . (A4) 

Here the variables {a{t) = a+{t) - 1, a{t),b{t) = b+{t) - 
l,b(t)} are classical quantities which correspond to the 
coherent-state eigenvalues, and the functional integral 



[4]) is performed subject to the conditions ai{oo) = 
ai(— oo) = bi{oo) = 6,;(— oo) = 0. The action S results 
from the Hamiltonian (IA2I) as 



S= dt 



- ^ (fli - aj){ai - a.j) 



<ij> 



1^ - ah{l + di)ai + PdiQi^l + bi)b 



i + IPidiai 



(A5) 

It is easily seen that the coupling induced by the num- 
ber 1 in the term abi{\ H- 0,^)0^ of the interaction part 
of (|A5p does not contribute to the calculation of the ex- 
pectation values (Greens functions) G'Ar({i, i}, {j, t'}) = 

(na=ia»»(*a)n^=iajf!(*^)>' especially to Pi{i,t) = 
Gi({i,t}, {0,0}), by perturbational series. Hence we ne- 
glect this destroying coupling in the following. Further- 
more, the variables bi and bi can be integrated out. E.g., 
performing the functional integral exp(— 5) leads 

to a factor ^ 5{^dtbi — (3diai{l+bi)) . These (5-conditions 
can be easily integrated (remember 6^(00) = 0) to 



6,(0 = cxp -/3 / dt' {h,ai){t') 



(A6) 



The remaining term in S which contains bi (t) yields then 



dtbi{t){diai){t) = aAi + — 6i(-oo) 

= ^Ai + 0{Al), (A7) 



where A, := ^_^dt{d.,ai){t). 

It remains to perform the average of the expectation 
values Gn over the disorder distribution of the pi . Note 
that the normalization factor of the path integral is de- 
fined so that / 2?[a, a, 6, 6] exp(— 5) = 1, and is purely 
numeric. Hence, their is no need for the replica trick. 
The average can be taken directly over exp(— 5). We use 
a Poissonian disorder distribution 



p{Pi = = exp(-p) 



(A8) 



of independent traps on each site i, characterized by the 
mean value p. We get 



exp 



fc=0 



7 / dtk{diai){t)^ 



exp |^p(^cxp(-7A) - 1 



exp 



PlAi 



^A? + 0(A^). 



(A9) 
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Using this expression, wc arrive at a reduced action for the calculation of the disorder averaged Greens functions 

dtl ^ a^dta^ + - ^ (a* - - aj) 

+ Y: [piA.. + ^^^A^ + 0{A^)] } . (AlO) 

i 

As long as a/3 — pj'^ > the third order terms 0{A^) 
becomes irrelevant in the RG sense. 

To obtain a proper field theoretic functional, it remains 
to transcribe the formulation from the lattice to the spa- 
tial continuum. Performing a naive continuum limit of 
S'rcd, several rescalings and a renaming of the variables 
and parameters, we finally arrive at the dynamical re- 
sponse functional p.4p . 



Appendix B: The renormalization of the longest SAW 



Here, we present some details of our diagrammatic calculation for the longest SAW. The general formula of the 
diagrammatic contributions related to the various SAWs we consider in this paper is 



/>oo 

/(q^r,wA)= / T[ds^D{q^,T,{si}exp[-iXwL{{si})] 
■Jo 



(Bl) 



where {s^} is the set of Schwinger-parameters of a given diagram. A is defined through Ao(A) = — iX^fLi ^ •= —iX- 
L{{si}) is a placeholder for the shortest length imin({si}), the longest length imax({si}) or the average length 
Lavi{si}), respectively, and w is a shorthand for the corresponding limit of Wr. Here, we will focus on Lmax({sj})- 
The average SAW and the corresponding multifractal moments are treated in the following appendix. 

First, let us focus on the 1-loop part of the calculation. For simplicity, we set q = in the following because we are 
not interested in reproducing the well-known field-renormalization Z. Moreover, we will neglect all contributions to 
diagrams proportional to r because we are not interested in reproducing the well-known Zr- At 1-loop order, there 
are only 2 conducting diagrams, namely diagrams A and B of Fig. ([IJ. These diagrams give 



5*0 = A - 2B = 



2 



2 

dsids2 



dsids2 



exp[-(si + S2)t] - 

exp[-(si + S2)t] 
[47r(si + sa)] 



cxp[—iXw max(si, S2)] — 2 exp[— iAwsi] 



d/2 



1 - 2^^5161(52 - Si) + 0(«/) 



Ge9^ -e/2 '''^'^ 



(B2) 



where w = W-q and where A is defined through A_o(A) = — «X]q=i 
(|6.8p . we get the w-part of the renormalized selfenergy 



r2li:'^ = ^f2 



(11) 



= iXwZ,, 



iX. Using the renormalization scheme 

2 s e/2-. 



iXwZwl 1 — 



{-h(i-i^)"ii(^r-°("')}' 



5 

12 



4£ 
2. e/2 



T 



It follows the renormalization factor Z^ to 1-loop order: 



Z„, = 1 + — + 0(u2) , 
4e 



(B3) 



(B4) 



Now, we turn to 2-loop order. First, we consider the diagrams C, D, E, F, G. These diagrams lead to the integral 



5i = C - 4D - E + 2F +4G = 



(47r)d 



Ql{{s,})cxp{-TJ2^S^) 



[{Sl + S2)(S3 -I- S4) -t- S5(S1 + S2 + S3 -I- S4)] 



d/2 ' 



(B5) 
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where we have used some invariance under permutations of indices to reduce the number of terms and where we have 
defined 

6i({si}) = exp[-zAw(si + S3)] ^0{si - S2 - S5)0{s3 - S4 - S5) - 26'(s3 - S4 - S5) - 0{si + S3 - S2 - S4) + 2| 

+ exp[-iAu;(si + S4 + S5)] |6'(.si + S5 - S2)0{si + 85- 83)9(31 +84,-82- S3) - 26'(s4 + S5 - S3) + l| . (B6) 
To simpUfy the integrations, we introduce new variables, 

Si — Xti , S2 = (1 — x)ti , 

S3 = yt2 , S4 = (1 - y)t2 , S5 = is . (B7) 
The integrations over x and ?/ are cumbersome but manageable and produce after expansion to linear order in w 

+ '-^S{ti - t2)9{t2 - ts) [3tl{t2 + ts) + 3t2tl +4 + 2tl] } . (B9) 
In the same manner we find for the second group of diagrams, H, I, J, K. L, 



W J [(s-^ +S5)(s2 +S4) +S2S4] 

with 

02({Si}) = exp[-iAti;s5] |i6'(s5 - si - S2 - S3)6'(s5 - si - S4 - S3) - 6'(s5 - si - S2 - S3) + i| 

+ exp[-iAw(si + S2 + S3)] |0(si + S2 + S3 - S5)^(s2 - S4) - 6'(si + S2 + S3 - 85) - 9{s2 - S4) + l| . 

(Bll) 

Here, we chose new integration variables 

8i = {l- x)yti, 83 = {I - x){l - y)ti , 

S2 = t2 , S4 = i3 , S5 = xti , (B12) 
and the integration over x and y yields 

Next we turn to the integration over the ^-variables. In part, these integrations can be done in an efficient and elegant 
manner by taking derivatives of the parameter integral 



6e 



a^"' (5 + c) + b^-' (a + c) + c^"' (a + 6)] - 3a6c| 

(B14) 



introduced by Breuer and Janssen [46| . The remaining parts can be tackled in the same spirit by introducing a second 
parameter integral: 

N{a, b, c) = -i-j / n dU ''J^^t^^lZ^^lZ^Qit, - t2)9{t2 - t3) . (B15) 

(4^) [tlt2+t2t3+t3hY'^ 
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In e-expansion, we obtain 
^(M,c) = - -(- 



a-'^b-^c] a^"^ + ( ^ - T In 3 ) 



35 1 

24 ^ 4 



i + iln3-iln2U^ 



12 4 



-§ + ^ln3-^ln2)a2&+(-l-^ln3+^ln2)a62 



31 3, „ 3, „\ , 

h -ln3+ - In 2 a^c 

16 8 4 ' 



i + ^ln3-^ln2)ac^ 



i - ^ln3+ |ln2^62c+ ( -\ + ^ In 3 - ^ In 2^ bc^ - ^abc 



.2 2 2 / V4 8 2 

Via differentiating the parameter integrals with respect to their parameters, we get the w-parts 



I2e 



■ iXw 



21 

ln2H hi3 



and 



S2 



12e 



■ iXw 



Collecting the 1- and 2-loop contributions, we find 



1(20 



rl 
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( 16e 48 
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8e 



r 15 



e \ T 



32e 



41 27 , 27 , 

n2H n3 

32 8 16 



49 35 , „ 69 , „ , 

ln2H ln3 u 

384 96 192 



9 ? 



£ \ T 



In 2 



(B16) 



(B17) 



(B18) 



(B19) 



for the w!-part of the renormalized self-energy to order 11} . This form makes evident that non-primitiv divergencies 
drop out and the e-poles are cancelled by choosing 



M / 15 3 

IH h \ 

4e V32e 128 



701n2-691n3 
192 



(B20) 



Appendix C: Renormalization of the multifractal moments 



Now, we present some details of our diagrammatic calculation for the multifractal moments. Our calculation here 
is based on Eq. (jBl[) with L({si}) specified to 



■yeB i 



(CI) 



where L^{{si}) = J2ie-y ^^'^ length of SAW 7 of the bundle B of all SAWs on the conducting part of the diagram, 
^(7) its probability, rrii the statistical weight of line i of the diagram. More generally, we are interested in all averaged 
moments of the weights, and hence we consider 



(C2) 



with a left general. For the 1- and 2-loop diagrams A to L shown in Fig. (1) the probabilities of the SAWs according 
to the kinetic rule are 

P({7})^^'={l/2,l/2}, p({7})(^)={l}, 

P({7})<''' = {1/4, 1/4, 1/4, 1/4} , p({7})(^) = Piil})^''^ = {1/2, 1/2} , 
M{7})'^'-M{7})'^'-{1}, 

PiM)^"^ = {1/2, 1/4, 1/4} , p({7})('^ = PiM)^-'^ = {1/2, 1/2} , 

P({7})<^'-M{7})'^^-{1}. (C3) 

Note that only the three SAWs on diagram H the kinetic rule lead to different probabilities than the static rule which 
yields pi{"/})stat ~ {1/3, 1/3, 1/3}. This fact is discussed in detail in the main text. The statistical weights of the 
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lines of the different diagrams follow from Eq. (|C1|) as 



1 

2 ' 



= 1 



4=1,2,3,4,5 — ™3,4,5 " "^1,2,3,4 



(F) _ (G) _ , 



'2,3,5 ~ '"2,3,5 



'2,4 



n,3,5 



2 ' 



'1,2,3,5 



2 ' 



"^1,3 = "^1,2,3 = "^5 - 1 



(C4) 

Using these weights as well as the symmetries of the diagrams, we obtain form Eq. (jC2p the following averaged 
moments of the weights: 



LtHA-2B,{s,}) = -[l-^){s, + s,) 
Li^J {C-AD-E + 2F + 4G, {s,}) = (l - ^) (^i 



32 + S3 + Si I 



4^)(i7-/-2J + 2i^ + L,{sJ)= fl 



2" 



.Si + S3 + S5 



1 - 



2-^)^5, 

1 \2, 



2" 



IS2 + S4) 



(C5) 



In the following, we replace the control parameter w by Va to emphasize the fact that the multifractal index a is 
kept general in our calculation. Using the substitutions si = a;t, S2 = (1 — x)t, the of the 1-loop self-energy that is 
linear in Va becomes 



'0,11 



-iXv. 



ds,ds2 ^-p;-(-^+;f] 4"j(A-2i.,{s.}) 



(si+S2)''/2 



e V 2" 

Using the renormalization scheme ()6.8p with the w renormalization replaced by 

Va ^ Va — ^ ^ ZaVa , 

we get 



(C6) 
(C7) 



= iX 



1 - 



1 \ U/M^\s/2 



1 



7 5 

2e " 12 



1 - 



2" / £ V r 



0(^3)} 



for the Wfj-part of the renormalized vertex-function to 1 -loop order. It follows 

1 \ u 



1+1 



2"/ e 



(C8) 



(C9) 



for the renormalization factors of the 

Now, we turn to 2-loop order. First, we consider the diagrams C, D, E, F, G. These diagrams lead to the integral 



Ii,v = -iXva 



9 



2{AttY 



-A- , {C~4D~E + 2F + 4G, {s,}) exp(-r J2^ Si) 

ll"^^*^^ — : — — \ — r — : — : — \ — \ — ^^7^ 



(CIO) 



i=l [(Sl + S2)(S3 + S4) + S5(S1 + S2 + S3 + Si)^ 

Switching to the integration variables defined in Eq. (|B7p and integrating over x and y, the integrals Ii^v can be 
expressed once again in terms of the mother- integral M . Taking the appropriate derivatives thereof, we obtain 

3 



h.v 



-iXvn 



9 



2{ATrY 
'2(1- 



1=1 
1 



ti<2exp[-T(fi +t2+ ta) 



[tlt2 + t2t3 + tsh 



d/2 



= —iXVa — 



(1 - ^) (*■+« + (2- |,)i3 
{('-^)f-(' 



o 

90+2 



)} 
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Proceeding similarly, we get for the second group of diagrams H, I, J, K, L 



-iXv, 



Lii {H-I-2J + 2K + i, {s J) exp(-T s. 



[(Sl + S3 + S5)(S2 + S4) + S2S4] 



d/2 



Using the integration variables defined in Eq. (jB12|) . we obtain 



"4(47r)'^ 



Y[dt, 



tl exp(-Tfti +t2 +^3)) 



d/2 



1 



l-^)*l + (l-i7)'(^2+t3) 



= —iXva — : — T 
As 



2" 
1 - 



1 



Collecting the 1- and 2-loop contributions, we find 



2^+1 )} 



(Cll) 



(C12) 



|(2/) 



(21) 

V 

1 

2" 



•M--(-^)[> + a-F-§)7]7(T) 

1 



iXVaS 1 



9 1 

2^2" 
I \ u 

1 



e/2 



15 



2c(+3 4Q + I 



9 9 



2"/ e 



9 47 

£ " 12 



11 

e 



2e2 
65\ 1 



0{u 



12 J 2" 



2 1\ 1 

e " 2/ 4" 



4e 



(C13) 



for the ti^-part of the renormalized self-energy to order u^. It is free of non-primitive divergencies as it should, and 
the £-poles are cancelled by the renormalization factors stated in Eq. (|7.10p . 
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